Exome sequencing identifies breast cancer susceptibility genes and defines the contribution of coding variants to breast cancer risk.

Linkage and candidate gene studies have identified several breast cancer susceptibility genes, but the overall contribution of coding variation to breast cancer is unclear. To evaluate the role of rare coding variants more comprehensively, we performed a meta-analysis across three large whole-exome sequencing datasets, containing 26,368 female cases and 217,673 female controls. Burden tests were performed for protein-truncating and rare missense variants in 15,616 and 18,601 genes, respectively. Associations between protein-truncating variants and breast cancer were identified for the following six genes at exome-wide significance (P < 2.5 × 10-6): the five known susceptibility genes ATM, BRCA1, BRCA2, CHEK2 and PALB2, together with MAP3K1. Associations were also observed for LZTR1, ATR and BARD1 with P < 1 × 10-4. Associations between predicted deleterious rare missense or protein-truncating variants and breast cancer were additionally identified for CDKN2A at exome-wide significance. The overall contribution of coding variants in genes beyond the previously known genes is estimated to be small.

A Comprehensive Analysis of 21 Actionable Pharmacogenes in the Spanish Population: From Genetic Characterisation to Clinical Impact.

The implementation of pharmacogenetics (PGx) is a main milestones of precision medicine nowadays in order to achieve safer and more effective therapies. Nevertheless, the implementation of PGx diagnostics is extremely slow and unequal worldwide, in part due to a lack of ethnic PGx information. We analysed genetic data from 3006 Spanish individuals obtained by different high-throughput (HT) techniques. Allele frequencies were determined in our population for the main 21 actionable PGx genes associated with therapeutical changes. We found that 98% of the Spanish population harbours at least one allele associated with a therapeutical change and, thus, there would be a need for a therapeutical change in a mean of 3.31 of the 64 associated drugs. We also identified 326 putative deleterious variants that were not previously related with PGx in 18 out of the 21 main PGx genes evaluated and a total of 7122 putative deleterious variants for the 1045 PGx genes described. Additionally, we performed a comparison of the main HT diagnostic techniques, revealing that after whole genome sequencing, genotyping with the PGx HT array is the most suitable solution for PGx diagnostics. Finally, all this information was integrated in the Collaborative Spanish Variant Server to be available to and updated by the scientific community.

A crowdsourcing database for the copy-number variation of the Spanish population.

BACKGROUND: Despite being a very common type of genetic variation, the distribution of copy-number variations (CNVs) in the population is still poorly understood. The knowledge of the genetic variability, especially at the level of the local population, is a critical factor for distinguishing pathogenic from non-pathogenic variation in the discovery of new disease variants. RESULTS: Here, we present the SPAnish Copy Number Alterations Collaborative Server (SPACNACS), which currently contains copy number variation profiles obtained from more than 400 genomes and exomes of unrelated Spanish individuals. By means of a collaborative crowdsourcing effort whole genome and whole exome sequencing data, produced by local genomic projects and for other purposes, is continuously collected. Once checked both, the Spanish ancestry and the lack of kinship with other individuals in the SPACNACS, the CNVs are inferred for these sequences and they are used to populate the database. A web interface allows querying the database with different filters that include ICD10 upper categories. This allows discarding samples from the disease under study and obtaining pseudo-control CNV profiles from the local population. We also show here additional studies on the local impact of CNVs in some phenotypes and on pharmacogenomic variants. SPACNACS can be accessed at: http://csvs.clinbioinfosspa.es/spacnacs/ . CONCLUSION: SPACNACS facilitates disease gene discovery by providing detailed information of the local variability of the population and exemplifies how to reuse genomic data produced for other purposes to build a local reference database.

A Large Case-Control Study Performed in Spanish Population Suggests That RECQL5 Is the Only RECQ Helicase Involved in Breast Cancer Susceptibility.

Around 50% of the familial breast cancer (BC) cases are estimated to be caused by germline variants in known low-, moderate-, and high-risk susceptibility genes, while the other half is of unknown genetic origin. In the present study, we wanted to evaluate the role of the RECQ helicases, some of which have been studied in the past as candidates, with unclear results about their role in the disease. Using next-generation sequencing (NGS) technology, we analyzed the whole coding sequence of BLM, RECQL1, RECQL4, RECQL5, and WRN in almost 2000 index cases from BC Spanish families that had previously tested negative for the known BC susceptibility genes (BRCAX) and compared the results with the controls extracted from gnomAD. Our results suggest that BLM, RECQL1, RECQL4, and WRN do not play a major role in BC susceptibility. However, in the combined analysis, joining the present results with those previously reported in a series of 1334 BC Spanish patients and controls, we found a statistically significant association between Loss of Function (LoF) variants in RECQL5 and BC risk, with an OR of 2.56 (p = 0.009; 95% CI, 1.18-4.98). Our findings support our previous work and places the RECQL5 gene as a new moderate-risk BC gene.

Genetic, environmental and life-style factors associated with longevity. Protocol and response of the LONGECYL Study / Genética, ambiente y estilo de vida asociados a longevidad. Protocolo y respuesta del estudio LONGECYL

Objective: To describe the objectives, the methodological approach, the response rate of the Genetic, Environmental and Life-style Factors Study in Castilla y León (Spain). Method: The Health Sentinel Network studied a sample of long-lived individuals aged 95 or more (LLI). The study included biological samples processed with the Global Screening Array v3.0 that contains a total of 730,059 markers. Written consent was obtained before the examination. Conclusions: The LLI contacted were 944, and 760 were completed studied. The 87.4% of LLI were born in Castile and Leon and only 1% were non-native of Spain. Severe cognitive impairment was declared in 8.1% of men and 19.2% of women. Genotyping was performed in 739 LLI, the 78.3% of the contacted sample. Family doctors and nurses achieve high participation in population-based studies. DNA samples were taken from 94% of fully studied LLI, and 100% of these samples where successfully genotyped. (AU)

Objetivo: Describir los objetivos, la metodología y la respuesta del estudio «Factores genéticos, ambientales y de estilo de vida asociados a la longevidad en Castilla y León» (España). Método: La Red Centinela Sanitaria estudió una muestra de personas longevas de 95 y más años de edad. Se pidió el consentimiento informado y se tomaron muestras biológicas para genotipificación, que fueron procesadas con el Global Screening Array v3.0 que contiene 730.059 marcadores. Conclusiones: Fueron contactados 944 longevos y 760 completaron el estudio. El 87,4% de los longevos habían nacido en Castilla y León, y solo el 1,5% no eran naturales. El 8,1% de los hombres y el 19,2% de las mujeres presentaban un deterioro cognitivo grave. Los médicos de familia y las enfermeras consiguen una alta participación y una buena calidad de la información en estudios de base poblacional. Se tomaron muestras de ADN en el 94% de los longevos y el 100% fueron genotipificadas con éxito. (AU)

Genetic, environmental and life-style factors associated with longevity. Protocol and response of the LONGECYL Study.

OBJECTIVE: To describe the objectives, the methodological approach, the response rate of the Genetic, Environmental and Life-style Factors Study in Castilla y León (Spain). METHOD: The Health Sentinel Network studied a sample of long-lived individuals aged 95 or more (LLI). The study included biological samples processed with the Global Screening Array v3.0 that contains a total of 730,059 markers. Written consent was obtained before the examination. CONCLUSIONS: The LLI contacted were 944, and 760 were completed studied. The 87.4% of LLI were born in Castile and Leon and only 1% were non-native of Spain. Severe cognitive impairment was declared in 8.1% of men and 19.2% of women. Genotyping was performed in 739 LLI, the 78.3% of the contacted sample. Family doctors and nurses achieve high participation in population-based studies. DNA samples were taken from 94% of fully studied LLI, and 100% of these samples where successfully genotyped.

POLRMT as a Novel Susceptibility Gene for Cardiotoxicity in Epirubicin Treatment of Breast Cancer Patients.

Anthracyclines are among the most used chemotherapeutic agents in breast cancer (BC). However their use is hampered by anthracycline-induced cardiotoxicity (AIC). The currently known clinical and genetic risk factors do not fully explain the observed inter-individual variability and only have a limited ability to predict which patients are more likely to develop this severe toxicity. To identify novel predictive genes, we conducted a two-stage genome-wide association study in epirubicin-treated BC patients. In the discovery phase, we genotyped over 700,000 single nucleotide variants in a cohort of 227 patients. The most interesting finding was rs62134260, located 4kb upstream of POLRMT (OR = 5.76, P = 2.23 × 10-5). We replicated this association in a validation cohort of 123 patients (P = 0.021). This variant regulates the expression of POLRMT, a gene that encodes a mitochondrial DNA-directed RNA polymerase, responsible for mitochondrial gene expression. Individuals harbouring the risk allele had a decreased expression of POLRMT in heart tissue that may cause an impaired capacity to maintain a healthy mitochondrial population in cardiomyocytes under stressful conditions, as is treatment with epirubicin. This finding suggests a novel molecular mechanism involved in the development of AIC and may improve our ability to predict patients who are at risk.

Breast Cancer Risk Genes - Association Analysis in More than 113,000 Women.

BACKGROUND: Genetic testing for breast cancer susceptibility is widely used, but for many genes, evidence of an association with breast cancer is weak, underlying risk estimates are imprecise, and reliable subtype-specific risk estimates are lacking. METHODS: We used a panel of 34 putative susceptibility genes to perform sequencing on samples from 60,466 women with breast cancer and 53,461 controls. In separate analyses for protein-truncating variants and rare missense variants in these genes, we estimated odds ratios for breast cancer overall and tumor subtypes. We evaluated missense-variant associations according to domain and classification of pathogenicity. RESULTS: Protein-truncating variants in 5 genes (ATM, BRCA1, BRCA2, CHEK2, and PALB2) were associated with a risk of breast cancer overall with a P value of less than 0.0001. Protein-truncating variants in 4 other genes (BARD1, RAD51C, RAD51D, and TP53) were associated with a risk of breast cancer overall with a P value of less than 0.05 and a Bayesian false-discovery probability of less than 0.05. For protein-truncating variants in 19 of the remaining 25 genes, the upper limit of the 95% confidence interval of the odds ratio for breast cancer overall was less than 2.0. For protein-truncating variants in ATM and CHEK2, odds ratios were higher for estrogen receptor (ER)-positive disease than for ER-negative disease; for protein-truncating variants in BARD1, BRCA1, BRCA2, PALB2, RAD51C, and RAD51D, odds ratios were higher for ER-negative disease than for ER-positive disease. Rare missense variants (in aggregate) in ATM, CHEK2, and TP53 were associated with a risk of breast cancer overall with a P value of less than 0.001. For BRCA1, BRCA2, and TP53, missense variants (in aggregate) that would be classified as pathogenic according to standard criteria were associated with a risk of breast cancer overall, with the risk being similar to that of protein-truncating variants. CONCLUSIONS: The results of this study define the genes that are most clinically useful for inclusion on panels for the prediction of breast cancer risk, as well as provide estimates of the risks associated with protein-truncating variants, to guide genetic counseling. (Funded by European Union Horizon 2020 programs and others.).

CSVS, a crowdsourcing database of the Spanish population genetic variability.

The knowledge of the genetic variability of the local population is of utmost importance in personalized medicine and has been revealed as a critical factor for the discovery of new disease variants. Here, we present the Collaborative Spanish Variability Server (CSVS), which currently contains more than 2000 genomes and exomes of unrelated Spanish individuals. This database has been generated in a collaborative crowdsourcing effort collecting sequencing data produced by local genomic projects and for other purposes. Sequences have been grouped by ICD10 upper categories. A web interface allows querying the database removing one or more ICD10 categories. In this way, aggregated counts of allele frequencies of the pseudo-control Spanish population can be obtained for diseases belonging to the category removed. Interestingly, in addition to pseudo-control studies, some population studies can be made, as, for example, prevalence of pharmacogenomic variants, etc. In addition, this genomic data has been used to define the first Spanish Genome Reference Panel (SGRP1.0) for imputation. This is the first local repository of variability entirely produced by a crowdsourcing effort and constitutes an example for future initiatives to characterize local variability worldwide. CSVS is also part of the GA4GH Beacon network. CSVS can be accessed at: http://csvs.babelomics.org/.

Regulatory CDH4 Genetic Variants Associate With Risk to Develop Capecitabine-Induced Hand-Foot Syndrome.

Capecitabine-induced hand-foot syndrome (CiHFS) is a common dermatological adverse reaction affecting around 30% of patients with capecitabine-treated cancer, and the main cause of dose reductions and chemotherapy delays. To identify novel genetic factors associated with CiHFS in patients with cancer, we carried out an extreme-phenotype genomewide association study in 166 patients with breast and colorectal capecitabine-treated cancer with replication in a second cohort of 85 patients. We discovered and replicated a cluster of four highly correlated single-nucleotide polymorphisms associated with susceptibility to CiHFS at 20q13.33 locus (top hit = rs6129058, hazard ratio = 2.40, 95% confidence interval = 1.78-3.20; P = 1.2 × 10-8 ). Using circular chromosome conformation capture sequencing, we identified a chromatin contact between the locus containing the risk alleles and the promoter of CDH4, located 90 kilobases away. The risk haplotype was associated with decreased levels of CDH4 mRNA and the protein it encodes, R-cadherin (RCAD), which mainly localizes in the granular layer of the epidermis. In human keratinocytes, CDH4 downregulation resulted in reduced expression of involucrin, a protein of the cornified envelope, an essential structure for skin barrier function. Immunohistochemical analyses revealed that skin from patients with severe CiHFS exhibited low levels of RCAD and involucrin before capecitabine treatment. Our results uncover a novel mechanism underlying individual genetic susceptibility to CiHFS with implications for clinically relevant risk prediction.

Association Between ABCB1 Genetic Variants and Persistent Chemotherapy-Induced Alopecia in Women With Breast Cancer.

Importance: Persistent chemotherapy-induced alopecia (pCIA) has been recently described in patients with breast cancer and in its most severe form occurs in up to 10% of these patients. Genetic risk factors associated with pCIA have not been adequately explored. Objective: To identify genetic variants associated with pCIA. Design, Setting, and Participants: In this genetic association study, 215 women with breast cancer treated with docetaxel-based chemotherapy with a follow-up of 1.5 to 10 years after the end of the treatment were recruited retrospectively through 3 hospital oncology units across Spain between 2005 and 2018. Severe pCIA was defined as lack of scalp hair recovery (Common Terminology Criteria for Adverse Events, version 3.0, grade 2) 18 months or more after the end of treatment. Patients with grade 2 pCIA were selected as cases, and those with no sign of residual alopecia 12 months after the end of docetaxel treatment were selected as controls. A genome-wide association study in a discovery phase was conducted, and logistic regression was used to identify variants associated with the risk to develop this adverse effect. The validity of the association was addressed through a replication phase. Exposures: Docetaxel-based chemotherapy. Main Outcomes and Measures: Genotypes of single-nucleotide variants associated with pCIA. Results: In total, 215 women with breast cancer (median age, 51.6 years; interquartile range, 44-60 years) were recruited (173 patients for the discovery phase and 42 patients for the replication phase). In the discovery phase, ABCB1 genetic variants were associated with risk to develop pCIA. In particular, single-nucleotide variation rs1202179, a regulatory variant located in an enhancer element that interacts with the ABCB1 promoter, was associated with the occurrence of pCIA. This finding was validated in the replication cohort (combined odds ratio, 4.05; 95% CI, 2.46-6.67; P = 3.946 × 10-8). This variant is associated with ABCB1 mRNA expression, and the risk allele was associated with decreased ABCB1 expression levels (P = 1.64 × 10-20). Conclusions and Relevance: This is the first study, to our knowledge, that identifies an association between a regulatory variant in the ABCB1 gene and the occurrence of pCIA in patients with breast cancer who were treated with docetaxel-based therapies. This finding suggests an important insight into the biological mechanisms underlying pCIA and opens the opportunity to explore personalized treatment of these patients.

BRAF V600E Detection in Liquid Biopsies from Pediatric Central Nervous System Tumors.

Pediatric Central Nervous System (CNS) tumors are the most fatal cancer diseases in childhood. Due to their localization and infiltrative nature, some tumor resections or biopsies are not feasible. In those cases, the use of minimally invasive methods as diagnostic, molecular marker detection, prognostic or monitoring therapies are emerging. The analysis of liquid biopsies which contain genetic information from the tumor has been much more widely explored in adults than in children. We compare the detection of BRAF V600E targetable mutation by digital-PCR from cell-free-DNA and EV-derived DNA (ctDNA) in serum, plasma and cerebrospinal fluid (CSF) isolated from a cohort of 29 CNS pediatric patients. Here we demonstrate that ctDNA isolated from serum and plasma could be successfully analyzed to obtain tumor genetic information which could be used to guide critical treatment decisions.

Corrigendum: miARma-Seq: a comprehensive tool for miRNA, mRNA and circRNA analysis.

This corrects the article DOI: 10.1038/srep25749.

The PNPLA3 Genetic Variant rs738409 Influences the Progression to Cirrhosis in HIV/Hepatitis C Virus Coinfected Patients.

Contradictory data about the impact of the rs738409 steatosis-related polymorphism within PNPLA3 gene on liver fibrosis progression in HIV/hepatitis C virus (HIV/HCV)-coinfected patients have been reported. Our objective was to test whether this, and other polymorphisms previously related to fatty liver disease in HIV infection linked to SAMM50 or LPPR4 genes, influence liver fibrosis progression in HIV/HCV-coinfected individuals. Three hundred and thirty two HIV/HCV-coinfected patients who consecutively attended four Spanish university hospitals from November 2011 to July 2013 were included. A liver stiffness cut-off of 14.6 kPa, as determined by transient elastography, was used to diagnose cirrhosis. Liver stiffness progression was studied in 171 individuals who had two available LS determinations without anti-HCV treatment between them. Moreover, 28 HIV/HCV-coinfected patients who underwent liver transplant, as well as 19 non-cirrhotic coinfected individuals used as controls, were included in an additional study. Only rs738409 was associated with cirrhosis: 45 (29.6%) of 152 G allele carriers versus 36 (20.0%) of 180 CC carriers showed cirrhosis (multivariate p = 0.018; adjusted odds ratio = 1.98; 95% confidence interval = 1.12-3.50). Also, 21 (30.4%) of 69 G allele carriers versus 16 (15.7%) of 102 CC patients showed significant liver stiffness progression (adjusted p-value = 0.015; adjusted odds ratio = 2.89; 95% confidence interval = 1.23-6.83). Finally, the proportion of rs738409 G allele carriers was significantly higher in transplanted individuals than in controls (p = 0.044, odds ratio = 3.43; 95% confidence interval = 1.01-11.70). Our results strongly suggest that the rs738409 polymorphism is associated with liver fibrosis progression in HIV/HCV-coinfected patients.

miARma-Seq: a comprehensive tool for miRNA, mRNA and circRNA analysis.

Large-scale RNAseq has substantially changed the transcriptomics field, as it enables an unprecedented amount of high resolution data to be acquired. However, the analysis of these data still poses a challenge to the research community. Many tools have been developed to overcome this problem, and to facilitate the study of miRNA expression profiles and those of their target genes. While a few of these enable both kinds of analysis to be performed, they also present certain limitations in terms of their requirements and/or the restrictions on data uploading. To avoid these restraints, we have developed a suite that offers the identification of miRNA, mRNA and circRNAs that can be applied to any sequenced organism. Additionally, it enables differential expression, miRNA-mRNA target prediction and/or functional analysis. The miARma-Seq pipeline is presented as a stand-alone tool that is both easy to install and flexible in terms of its use, and that brings together well-established software in a single bundle. Our suite can analyze a large number of samples due to its multithread design. By testing miARma-Seq in validated datasets, we demonstrate here the benefits that can be gained from this tool by making it readily accessible to the research community.

Hepatitis C virus reinfection after sustained virological response in HIV-infected patients with chronic hepatitis C.

OBJECTIVES: To assess the incidence of hepatitis C virus (HCV) reinfections after therapy-induced clearance in HIV-coinfected patients with prior chronic hepatitis C. METHODS: Eighty-four HIV-infected subjects, who had previously achieved sustained virological response (SVR) after being treated of chronic hepatitis C, were analyzed. In all of them, at least yearly HCV RNA determinations were carried out during a median (range) of 34 (12-146) months. RESULTS: Seventy-two (86%) subjects had been people who inject drugs (PWID), of whom 11 (15%) continued to use snorted or injected drugs during the follow-up. Four (4.76%) patients showed HCV reinfection (incidence 1.21 [95% confidence interval: 0.3-3.09] cases per 100 person-years). These patients maintained risk factors for HCV infection. In three cases, HCV genotype switched. Phylogenetic analysis of the remaining case suggested reinfection from his sexual partner. CONCLUSION: The incidence of HCV reinfection in the overall population of HIV-coinfected patients who achieved SVR after being treated against chronic hepatitis C is low. A low frequency of risk behavior is the main factor accounting for this modest rate of reinfection. The possibility of reinfection should not be considered a reason against treatment of HCV infection with direct acting antivirals in PWID.

Reassessment of genotype 1 hepatitis C virus subtype misclassification by LiPA 2.0: implications for direct-acting antiviral treatment.

The accuracy of LiPA 2.0 for hepatitis C virus 1 (HCV-1) subtype classification was analyzed. LiPA 2.0 genotype results from 101 HCV-1-infected patients were compared to genotype findings determined by direct core sequencing. Eleven (11%) samples were misclassified. Given the influence of the HCV-1-subtype in the anti-HCV therapy response, an alternative classification method is warranted.

Contributions of PHOX2B in the pathogenesis of Hirschsprung disease.

Hirschsprung disease (HSCR) is a congenital malformation of the hindgut resulting from a disruption of neural crest cell migration during embryonic development. It has a complex genetic aetiology with several genes involved in its pathogenesis. PHOX2B plays a key function in the development of neural crest derivatives, and heterozygous mutations cause a complex dysautonomia associating HSCR, Congenital Central Hypoventilation Syndrome (CCHS) and neuroblastoma (NB) in various combinations. In order to determine the role of PHOX2B in isolated HSCR, we performed a mutational screening in a cohort of 207 Spanish HSCR patients. Our most relevant finding has been the identification of a de novo and novel deletion (c.393_410del18) in a patient with HSCR. Results of in silico and functional assays support its pathogenic effect related to HSCR. Therefore our results support that PHOX2B loss-of-function is a rare cause of HSCR phenotype.

Four new loci associations discovered by pathway-based and network analyses of the genome-wide variability profile of Hirschsprung's disease.

Finding gene associations in rare diseases is frequently hampered by the reduced numbers of patients accessible. Conventional gene-based association tests rely on the availability of large cohorts, which constitutes a serious limitation for its application in this scenario. To overcome this problem we have used here a combined strategy in which a pathway-based analysis (PBA) has been initially conducted to prioritize candidate genes in a Spanish cohort of 53 trios of short-segment Hirschsprung's disease. Candidate genes have been further validated in an independent population of 106 trios. The study revealed a strong association of 11 gene ontology (GO) modules related to signal transduction and its regulation, enteric nervous system (ENS) formation and other HSCR-related processes. Among the preselected candidates, a total of 4 loci, RASGEF1A, IQGAP2, DLC1 and CHRNA7, related to signal transduction and migration processes, were found to be significantly associated to HSCR. Network analysis also confirms their involvement in the network of already known disease genes. This approach, based on the study of functionally-related gene sets, requires of lower sample sizes and opens new opportunities for the study of rare diseases.